Universal thermodynamics of a two-dimensional Bose gas 
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Using renormalization-group arguments we show that the low-temperature thermodynamics of 
a three- or two-dimensional dilute Bose gas is fully determined by a universal scaling function 
Td(nl ksT, g(T)) once the mass m and the s-wave scattering length ad of the bosons are known (d 
is the space dimension). Here \i and T denote the chemical potential and temperature of the gas, 
and the temperature-dependent dimensionless interaction constant g(T) is a function of ma^/c^T/fr 2 . 
We compute the scaling function Ti using a nonperturbative renormalization-group approach and 
find that both the fi/ksT and g(T) dependencies are in very good agreement with recent exper- 
imental data obtained for a quasi-two-dimensional Bose gas with or without optical lattice. We 
also show that the nonperturbative renormalization-group estimate of the Berezinskii-Kosterlitz- 
Thouless transition temperature compares well with the result obtained from a quantum Monte 
Carlo simulation of an effective classical field theory. 

PACS numbers: 05.30.Rt,05.30.Jp,67.85.-d,03.75.Hh 



I. INTRODUCTION 

The physical properties of a system at thermal equilib- 
rium are determined by an equation of state. For a fluid 
of particles in the grand canonical ensemble, the equa- 
tion of state relates a thermodynamic quantity such as 
pressure, density or entropy to temperature and chemi- 
cal potential. It can take a complicated expression when 
the particles interact via a two-body potential V(ri — r 2 ) 
which has no simple expression as is usually the case 
in real systems. Quite remarkably however, the low- 
temperature equation of state of a d-dimensional dilute 
gas is universal, in the sense that it depends only on 
a small number of parameters, such as the mass m of 
the particles and the s-wave scattering length a^, and is 
otherwise insensitive to the details of the two-body po- 
tential V(ri — r2). A well-know example of universality 
is given by a three-dimensional dilute Bose gas at zero 
temperature, the pressure of which is given by the mean- 
field result P{fj) = /j J 2 m/(87rash 2 ) to leading order in the 
small parameter ma 2 fi/h 2 . The first quantum correction, 
known as the Lee-Huang- Yang correction, is also entirely 
determined by m and as (besides the chemical potential 
/i) [ ]. (In the following we set ks = h = 1.) 

The universality property of the equation of state of 
a dilute Bose gas can be understood from the point of 
view of the theory of phase transitions [2-5] . By varying 
the chemical potential from negative to positive values at 
zero temperature, one induces a quantum phase transi- 
tion between a state with vanishing pressure and no parti- 
cles (vacuum) and a superfluid state with a nonzero pres- 
sure. This identifies the point /i = T = as a quantum 
critical point (QCP). Above two dimensions (the upper 
critical dimension of the T = quantum phase transi- 
tion), i.e. for d > 2, the boson-boson interaction is irrel- 
evant and the critical behavior at the transition is mean- 
field like with a correlation-length exponent v = 1/2 and 
a dynamical exponent z — 2. However, the boson-boson 



interaction cannot be completely ignored and enters the 
equation of state [ ]. In the critical regime near the QCP, 
defined by ml 2 \± <C 1 and ml 2 T <C 1 with I the "natural" 
low-energy length scale [ , ] , the pressure takes the form 



%r) = g) d/2 T^j,(^(r) 



(i) 



where T& is a universal scaling function characteris- 
tic of the (i-dimensional dilute Bose gas universality 
class. The temperature-dependent dimensionless inter- 
action constant g(T) is a known function of ma^T, so 
that P([i, T) can also be written in terms of a universal 
function of \ijT and ma 2 d T. In two dimensions and in the 
weak- interaction limit, g(T) = g is approximately tem- 
perature independent and the universal scaling function 
^(fi/T, g) depends on ji/T with the interaction strength 
g as a parameter; the equation of state then exhibits an 
approximate scale invariance (with no characteristic en- 
ergy scales other than \i and T) [9, 10]. Equation (1) 
also holds in a one-dimensional Bose gas (i.e. below the 
upper critical dimension of the T — vacuum-superfluid 
transition) but with the universal function T\ depending 
only on ji/T. 

While Eq. (1) follows from general renormalization- 
group (RG) arguments (see Sec. Ill), the theoretical de- 
termination of the universal scaling function J^(x, y) re- 
quires an explicit computation of the pressure P(/i, T). A 
perturbative calculation order by order in g(T) is possible 
only for d > 2 (it nevertheless breaks down in the critical 
regime of the thermal phase transition between the nor- 
mal and the superfluid phase, which is controlled by the 
Wilson-Fisher fixed point of the classical 0(2) model). 
In two dimensions, perturbative theory is plagued with 
infrared divergences at finite temperatures, thus mak- 
ing the determination of T^ difficult, in particular in 
the quantum critical regime \\i\ <C T. The Berezinskii- 
Kosterlitz-Thouless (BKT) transition [ ] and the low- 
temperature phase with quasi-long-range order are also 
beyond a mere perturbative treatment [12]. 



The advantage of the point of view based on phase 
transitions is two-fold. Firstly it gives a straightforward 
explanation of universality in a dilute Bose gas. Sec- 
ondly it shows that the universal equation of state (1) 
holds not only for a dilute Bose gas but for any system 
near a quantum phase transition belonging to the same 
universality class. For instance, a Bose gas in an opti- 
cal lattice near the vacuum-superfluid transition exhibits 
the same thermodynamics as a dilute Bose gas, provided 
that m and a^ are understood as the effective mass and 
scattering length of the bosons moving in the lattice. The 
thermodynamics of a Bose gas near the superfluid-Mott- 
insulator transition is also described by the equation of 
state (1), since this quantum phase transition (when it is 
induced by a density change) belongs to the dilute Bose 
gas universality class [2, 3]. In this manuscript we focus 
on the vacuum-superfluid transition in a two-dimensional 
Bose gas. 

On the experimental side, cold atomic gases provide us 
with highly controlled and tunable systems where uni- 
versal thermodynamics can be experimentally demon- 
strated. Altough cold gases are inhomogeneous and of 
finite size due to the harmonic confining potential, us- 
ing a local-density approximation it is possible to deduce 
the equation of state P([i,T) of the infinite homogeneous 
gas (with uniform density) [13, 14]. A number of experi- 
ments on weakly-interacting two-dimensional Bose gases 
have been reported [ - ], and the scale invariance of 
the equation of state P([i,T) has been observed [ 7-19]. 
More recently, the equation of state of a Bose gas in 
an optical lattice has been measured near the vacuum- 
superfluid transition in a regime where the interaction 
constant is not weak [' ]. These experiments allow us 
to determine both the ji/T and g(T) dependence of the 
universal scaling function Ti in various limits and will be 
thoroughly discussed in the manuscript. 

The outline of the paper is as follows. In Sec. II, we 
introduce and motivate the low-energy effective Hamilto- 
nians which enable to derive the universal thermodynam- 
ics of three- and two-dimensional dilute Bose gases. Sec- 
tion III is devoted to a discussion of the thermodynamics 
of a dilute Bose gas using the language and concepts fa- 
miliar from the theory of phase transitions. A detailed 
derivation of Eq. (1) is given. In Sec. IV, we discuss 
the universal scaling function Ti obtained from a non- 
pert urbative renormalizat ion- group (NPRG) approach. 
These theoretical results are compared with the exper- 
imental data of Refs. [18-20] in Sec. V. In particular, we 
make quantitative comparisons between the experimen- 
tal data obtained with a Bose gas in an optical lattice [2 } 
and theoretical results obtained in the framework of the 
Bose-Hubbard model. Finally, in Sec. VI, we discuss the 
NPRG prediction for the BKT transition temperature 
and compare it with the estimate deduced from a quan- 
tum Monte Carlo simulation of an effective classical field 
theory. 



II. MODEL HAMILTONIANS 



A. Three-dimensional Bose gas 

The interaction between ultracold atoms is governed 
by a potential V( r i — r^) which is repulsive at short dis- 
tances and determined by the van der Waals attraction 
— Cq\ti — T2\~ 6 at long distances [ , ]. The latter de- 
fines the microscopic length scale / v dw ~ (mCe) 1 ^ (m 
denotes the atomic mass). For length scales larger than 
/ v dw and energies smaller than l/m^ dW , collisions be- 
tween atoms occur only in the s-wave channel and the 
scattering amplitude is well approximated by 



/3D(q) = 



a 3 



1 + z|q|a 3 ' 



(2) 



where the three-dimensional s-wave scattering length a 3 
is typically of the order of Z v dW- I n this low-energy 
regime, the ultracold gas can be described by the effective 
Hamiltonian 



H 



Jd d r^(r)(- 



2m 



[i ip(r) 



IM 



^T(r)^T(r)^(r)^(r) 



(3) 



3 for 



with an ultraviolet momentum cutoff A ~ ^ v dw (^ 
a three-dimensional gas). Here ^(^(r) is a bosonic op- 
erator, jjl the chemical potential and g the "microscopic" 
interaction constant. The scattering amplitude obtained 
from (3) takes the form (2) with a scattering length 



a 3 



mg 



47r + ^mgA 



(4) 



which is a function of g and A. The low-energy effec- 
tive description is valid only for momentum scales much 
smaller than A, which requires both temperature and 
density to be small enough: T <C A 2 /2m and D <C A 3 . 



B. Two-dimensional Bose gas 

A quasi-two-dimensional gas can be created by sub- 
jecting a three-dimensional gas to a confining harmonic 
potential along one direction. The scattering amplitude 
then vanishes in the low-energy limit q — >• 0, 



/2D(q) = -- 



2tt 



In 



(J^) 



c 



(5) 



(C is the Euler constant), as in a strictly two-dimensional 
system [22, 23]. The effective two-dimensional s-wave 
scattering length a2 is a function of the thickness l z of 
the gas in the confining direction, as well as the s-wave 
scattering length and microscopic interaction strength of 
the three-dimensional (unconfined) Bose gas. At suffi- 
ciently low temperatures, when T is much smaller than 



the uo z = l/ml 2 z , only the lowest level of the confining 
potential is populated and the gas behaves as a two- 
dimensional system. The quasi-two-dimensional gas can 
be described by the effective Hamiltonian (3) with d = 2 
and a "microscopic" interaction constant [22, 23] 



ml z 



(6) 



which reproduces the scattering amplitude (5) with the 
scattering length 



CL2 



A 



exp 



2tt 
mg 



C 



(7) 



Here A ~ l~ x is an ultraviolet momentum cutoff below 
which the two-dimensional description holds. In addition 
to the condition T <C uj z ~ A 2 /2m, we must require 
the density to satisfy D « A 2 . The typical energy per 
particle gD is then much smaller than uj z as it should be 
for the two-dimensional description to be justified. Note 
that in all experiments realized so far, the dimensionless 
interaction constant g = 2mg is small. 



C. Bose gas in an optical lattice 

Bosons in an optical lattice are described by the Bose- 
Hubbard model [3, 241, 



H 



<r,r') 



(ftjr 



h.C.)+£ 



fin r + —n r (n r - 1) 



(8) 

where t is the hopping amplitude between nearest- 
neighbor sites (r, r'), and U the onsite interaction, fy 
is an annihilation (creation) operator for a boson at site 

r of the lattice and n r = V4^ r - ^ n e ^ ec ^ ve single- 
band description is valid only if the optical potential 
is strong enough and at sufficiently low temperatures. 
For a (i-dimensional hypercubic lattice, the dispersion 
of the free bosons is given by the Fourier transform 
t q = — 2t ^2 i cos(qil) of the intersite hopping matrix (I 
denotes the lattice spacing). It is convenient to use a 
shifted dispersion law 



d 



2td-2t) cosfeZ) 



(9) 



i=l 



and behaves as e c 



tl 2 q 2 for 



which vanishes for q 
|q|«r 1 . 

If the density D is low enough (Dl d < 1), the ground 
state is always a superfluid and we do not have to worry 
about the physics of the Mott transition [3]. Further- 
more, at low temperatures T< t, the lattice does not 
matter and one can take the continuum limit where the 
Hamiltonian takes the form (3) with an effective mass 
m = l/2tZ 2 , an interaction constant g = Ul d , and a chem- 
ical potential \i + 2dt. The ultraviolet momentum cutoff 



A is of the order of the inverse lattice spacing Z _1 ; the 
conditions T <C t and Dl d <C 1 then become T <C A 2 /2m 
and D <C A d . 

Thus, in the low-energy limit, a Bose gas in an optical 
lattice behaves similarly to a homogeneous Bose gas with 
an effective mass m and an effective interaction constant 
g. To ensure that the effective continuum model repro- 
duces the same low-energy physics as the lattice model, 
we must choose the cutoff A so that it yields the same 
scattering length. In the two-dimensional case, we re- 
quire Eq. (7) to reproduce the scattering length of the 
two-dimensional Bose-Hubbard model [251, 



a 2 



2v^ 
which gives A = Ay/2 /I. 



exp 



Ant 
IT 



C 



(10) 



III. UNIVERSAL THERMODYNAMICS 

In this section, we discuss the thermodynamics of a 
d-dimensional dilute Bose gas from the point of view 
of phase transitions, starting from the Hamiltonian (3). 
This description provides us with a natural explanation 
of universality as well as a simple derivation of Eq. (1). 
Altough we will mainly focus on two-dimensional sys- 
tems in the following sections, for generality we consider 
an arbitrary dimension d>2. 



A. Vacuum-superfluid transition 

Let us first consider the vacuum-superfluid quantum 
phase transition induced by a change of chemical poten- 
tial at zero temperature. For d larger than the upper 
critical dimension d+ = 2, boson-boson interactions are 
irrelevant (in the RG sense) and the critical behavior is 
described by non-interacting bosons. At the QCP \i = 0, 
the ground state is the vacuum, and the single-particle 
Green function is given by 



G(q,i(j) 



2m 



(11) 



with uj a (bosonic) Matsubara frequency. This result is 
exact and holds for any value of the (bare) interaction 
constant g p ]. We deduce the dynamic exponent z = 2 
while the anomalous dimension r] vanishes. Similarly, for 
/i < 0, we find G(q, iuo) -1 = iu + /i — q 2 /2ra and the 
critical exponent associated with the correlation length 
£ = |2m/i| -1 / 2 takes the value v = 1/2. The value of 
the renormalized interaction gn at the QCP is given by 
the T matrix in vacuum (using again the fact that the 
ground state is the vacuum). In the low-energy limit, it 
takes the value Airas/m in three dimensions, but vanishes 
logarithmically in two dimensions (see Eq. (14) below). 



The same analysis holds for bosons moving in a lattice. 
In the vacuum, the single-particle Green function is given 

by 



G(q, iuo) = (ioj + \i + 2dt — e q ) 



(12) 



where e q is the dispersion of the free bosons [Eq. (9)]. 
The T = QCP between the vacuum and the superfluid 
phase is now located at \i c = —2dt and the elementary 
excitations have an effective mass m = l/2tl 2 . As in the 
continuum model, the renormalized value qr of the inter- 
action (i.e. the T matrix in vacuum) can be expressed in 
terms in of the scattering length a^ of the bosons moving 
in the lattice [27]. 



B. RG approach 

The preceding results can be formulated in the lan- 
guage of the RG. In the Wilson formulation, a RG trans- 
formation consists in integrating out "fast" modes with 
momenta between A and A/5 (s > 1), and rescaling fields, 
momenta and frequencies in order to restore the original 
value of the cutoff A. This yields an effective Hamiltonian 
for the "slow" modes with a renormalized interaction con- 
stant g(s) [ ]. At the QCP /i = T = 0, there is no 
renormalization of the quadratic part of the Hamiltonian, 
in agreement with the fact that Eq. (11) is exact. The 
dimensionless interaction constant g(s) = 2mA d ~ 2 g(s) 
satisfies the RG equation 



ds 



d)g(s) 



Kd~, ,-2 



(13) 



where K d = [2 d ~ 1 7r d/2 T(d/2)]~ 1 . Above the upper crit- 
ical dimension d+ = 2, g(s) vanishes for s — >• 00 and 
the only fixed point of Eq. (13) is the Gaussian fixed 
point g = 0, which therefore governs the quantum phase 
transition between the vacuum and the superfluid phase. 
From (13), we obtain 



9(s) 



87rAa 3 



47T 



ln(4^)+C 



if d = 3, 

if d = 2, 



(14) 



where the result for d = 3 holds for s ^> 1. 

There are two relevant perturbations about the Gaus- 
sian fixed point \i = T = g = 0: the chemical poten- 
tial and the temperature. In a RG transformation, they 
transform as 



f(s) 



z f, £l(s) 



8 1/v fi, 



(15) 



near the QCP {i.e. when |/2(s)| < 1 and f(s) < 1). We 
have introduced the dimensionless variables [29] 



2mT 



fi 



2m [i 



(16) 




FIG. 1. (Color online) Phase diagram of the dilute Bose gas 
(d > 2). The dashed lines are defined by |/x| ~ T and the 
solid one corresponds to the superfluid transition (of BKT 
type when d = 2). The shaded area corresponds to the high- 
energy region |/x|,T > A 2 /2m where the thermodynamics is 
not universal. The value of the ultraviolet momentum cutoff A 
is discussed in Sec. II for a three- and a quasi-two-dimensional 
Bose gas. 



Note that in the low-temperature regime where this 
analysis based on the effective Hamiltonian (3) is valid, 
|/2| C 1 and T <C 1 (see Sec. II). When \i and T are 
nonzero, the RG equation for g(s) is well approximated 
by (13) or (14) as long as \£l(s)\ < 1 and f(s) < 1. We 
can obtain a rough sketch of the phase diagram by not- 
ing that the low-energy behavior of the system depends 
on which of the conditions |/2(s)| ~ 1 and T(s) ~ 1 is 
reached first. This yields two crossover lines defined by 
|/2 1 ~ T, i.e. |/i| ~ T using z = 1/v = 2, in agree- 
ment with the generic phase diagram of a system near a 
quantum critical point (see Fig. 1) [ ]. For \i < and 
\fi\ ^> T, the system behaves as a dilute classical gas and 
we expect a classical Boltzmann description to apply (see 
Sec. IV). The condition |/i| <C T defines the quantum crit- 
ical regime where the physics is controlled by the QCP 
[i — T — and its thermal excitations [ ] . 



C. Universal thermodynamics 

Let us now consider the dimensionless pressure [29] 

2m 



P(P>,T,g) 



k d +> 



:P(»,T), 



(17) 



expressed in terms of the dimensionless variables T, ft 
and g (note that P has no explicit dependence on the 
ultraviolet cutoff A). In a RG transformation, P trans- 
forms as 



P(j2,f,g) = s- 



*P{s l ' v il,s z f,~g{s)), 



(18) 



provided that T(s) <C 1 and \p>(s)\ <C 1. Equation (18) 
holds for the full pressure since the vanishing of P when 
fi < and T = implies that P has no regular part at 
the transition. Only the two-body interaction constant g 
is taken into account. Higher-order interactions (which 



are inevitably generated in the RG procedure), such as 
the three-body term, are not considered here since they 
are irrelevant and give rise to subleading contributions to 
the pressure [30, 31] 

Setting s = f' 1 ^ in Eq. (18), we obtain 



Pfaf,g) 



fd/z+l 

(4^) rf / 2 



T d 



fl/v; 



;9(f) 



(19) 



where we use the notation g(T) for g(s = T -1 / 2 ). Going 
back to dimensionful variables and setting z = 1/v = 
2, we finally obtain Eq. (1) where the energy-dependent 
interaction constant 



87Ty / zmap 



~9(e) = { 



47T 



In 



( \ ^2ma\ 



if d- 
if d- 



3, 
2, 



(20) 



C 



is obtained from (14) with s = e -1 / 2 and e = 2me/A 2 . 
We stress that Td is a universal scaling function charac- 
teristic of the d-dimensional dilute Bose gas universality 
class (the factor l/(47r) d / 2 in (19) is introduced for later 
convenience); it is independent of microscopic parame- 
ters such as the mass m of the bosons or the scattering 
length ad which depend on the system considered. 

Using Eq. (1), we can write any thermodynamic quan- 
tity in a scaling form. For instance, the density D = 
dP/dji and the entropy per unit volume s = dP/dT read 



fmT\ 



d/2 



?, 



(1.0) 



(£,scn) 



(21) 



and 



s(fi,T) 



~2^J 



/i 77(1.0) 

Jl^d 






where T 



(i,o) 



(x,y) 



d x T d {x,y), T^ 0,1 \x,y) 

d-2 



(22) 



g(T) 



dyTd(x, y), and g'(T) = dg/dT. Since Tg'(T) 

for d > 2 and Tg'(T) = g(T) 2 /8Tr for d = 2, s(^,T) is a 

function of [ijT and g(T) (up to the factor (mT) d / 2 ). 

One often introduces the so-called phase-space pres- 
sure and phase-space density, 



V(fi,T) = P(n,T) 
T>(ji,T) 



A dB 



T 
D(^T)\ d dB 






(23) 



T\ 



where AdB = \/2ii/mT is the thermal de Broglie wave- 
length. V and V provides a direct measure of the scaling 



function Td and its derivative T ( 



(1,0) 



One can also con- 



sider the entropy per particle S = D 1 dP/dT, 

T d 



SQi,T) 



El 
t 



i" 



-F 



(1,0) 



T-(O.l) 

Tg'(T)- d 



r 



(1,0) ; 



(24) 



where we use the shorthand notation Td = 
Fd{n/T,g{T)), etc. Note that S({i,T) is also a 
universal function of ji/T and g(T) (see the remark 
about Tg'(T) following Eq. (22)). 

At zero temperature, the scaling function Td can be 
computed in perturbation theory. The one-loop correc- 
tion to the mean-field result gives 



P(/i,0) = O(/i) 



87r<23 



1 



64 
15tt 



ma\\± 



(25) 



in three dimensions (the one-loop correction is known as 
the Lee-Huang- Yang correction [ ]), and 



F(/i,0) = -O(/i) 



47T 



In 



ma\ii 



C 



(26) 



in two dimensions, where denotes the step func- 
tion [33-35]. These results can be cast in the scaling 
form 



p<,,o = Ef 



d/2+1/ 



/i 



,#0) 



(27) 



which is equivalent to Eq. (1) but more appropriate to 
the zero-temperature limit [36]. 

At finite temperature, the determination of the scal- 
ing function Ti (or Q^) is difficult in two dimensions, in 
particular in the quantum critical regime |/i| <C T. In 
the following section, we discuss the scaling function T2 
obtained from the NPRG approach. 



IV. SCALING FUNCTION J 7 OF A 
TWO-DIMENSIONAL BOSE GAS 

The NPRG approach has recently been used to un- 
derstand the physics of a Bose gas beyond the Bogoli- 
ubov approximation [ 7-42], but the computation of the 
scaling function Td has not been carried out except for 
the Lee-Huang-correction in a zero-temperature three- 
dimensional Bose gas [40]. Here we discuss the NPRG 
results for the scaling function T = T2 which determines 
the thermodynamics of a two-dimensional Bose gas. We 
use both the standard version of the NPRG as well as 
its lattice version [27, 43] to directly study the Bose- 
Hubbard model. The NPRG approach is briefly reviewed 
in Appendix A. Our results are based on the numerical 
solution of the NPRG equations as well as analytical re- 
sults in some limits, in particular for \i — (see Ap- 
pendix B). 



A. T(x, y) vs x (y fixed) 

We first discuss the x dependence of T(x : y) for fixed 
y. Figure 2 shows the phase-space pressure V = T 
[Eqs. (23)] as a function of /i/T for g(T) = 0.22 and 
g(T) = 5. We can verify that the scaling form (1) 




FIG. 2. (Color online) Phase-space pressure V(/jl,T) vs 
fi/T at fixed g(T). The upper symbols are obtained for 
g(T) = 0.22 with T = 0.1A 2 /2ra and 2m = 1. (Red) 
squares: (T, g, ra, A), (black) dots: (T/2,g, ra, A), (orange) 
diamonds: (T, 2g,m/2, A). The lower symbols are obtained 
for </(T) = 5. (Blue) triangles: (T,g,m,A), (green) trian- 
gles: (T/2,g,m,A/y/2), (purple) triangles: (2T,2g,m/2,A). 
The solid lines are guides to the eyes. (/x/T)bkt at the BKT 
transition is given in table I. 



9(T) 




0.1 0.22 0.5 


1 


(/x/T)bkt 




0.08 0.15 0.29 0.51 


(/x/T)bkt (Ref. [ 


]) 


0.09 0.17 0.31 


0.51 



TABLE I. (/x/T)bkt at the BKT transition for various various 
values of g(T) < 1 as obtained from the NPRG (Sec. VI) and 
Monte Carlo simulations [10]. 



V 50 




holds by computing V for various sets of parameters 
(T, g,ra, A). For the case g(T) = 0.22, we choose 
the value g = 0.22 for the bare interaction constant, 
so that the system is in the weak-coupling limit and 
g(T) ~ 0.22 nearly temperature independent (see the 
discussion below). We find that the three sets of pa- 
rameters (T, g, ra, A), (T/2,#,ra, A) and (T, 2g,m/2, A) 
(with T = 0.1A 2 /2ra and 2m = 1) yield the same 
results for the phase-space pressure V in agreement 
with the expected scale invariance at weak coupling: 
r(n/T,g(T)) = F(lJL/T,g). In the case g(T) = 5, 
the results obtained for the three sets of parameters 
(T,#,ra,A), (T/2,g,m,A/y/2) and (2T, 2g,m/2, A) also 
collapse on a single curve corresponding to the scaling 
function F(x,y) (with y = g(T) fixed). In this case one 
must change simultaneously at least two parameters to 
keep g(T) unchanged. In table I we indicate the value 
of (fi/T) bkt at the BKT transition for various values of 
g(T) < las obtained from the NPRG (Sec. VI) and 
Monte Carlo simulations [ ]. Note that neither our 
method nor the Monte Carlo simulations gives a reli- 
able estimate of (/i/T) bkt m the strong-coupling limit 

g(T)>i. 

Figure 3 shows the phase-space pressure V = J 7 , the 



FIG. 3. (Color online) Phase-space pressure V(fi,T), phase- 
space density T)(/j,,T) and entropy per particle S(/jl,T) vs 
/i/T. From top to bottom, the curves correspond to g(T) = 
0.1,0.5,1,3,5 for V, g(T) = 0.5,1,3,5 for T>, and g(T) = 
0.1,0.5,1 for S. The dots show the limiting behavior (30) 
valid for \i ^> T and the diamonds the classical gas result (28). 
[For numerical reasons, it is difficult to compute the entropy 
when g(T) > 1.] 



phase-space density V = J 7 ^ 1 ' ) and the entropy per par- 
ticle S(n,T) as a function of ji/T for various values of 
g(T) [Eqs. (23-24)]. At large and negative chemical po- 
tential (|/i|/T ^> 1), we find that the system behaves as 
a classical dilute gas, 



which implies 

T[x,y) -- 



n^T) = 


P(/i,T) 


= e 


5(/x,T) = 2- 


~T' 


= e x for 


x < 


and 



X > 1. 



(28) 



(29) 



In the opposite limit of a large positive chemical potential 
(///T ^> 1), the pressure can be approximated by its zero- 



temperature limit. Using the expression (26), we obtain 
2 / 2tt 1 . 



T>{n,T) 



a ( 4ir 1 , „ 
= I — + -ln2 



In 2 



T \g(ji) 2 



(30) 



S(ji,T) = 0, 



i.e. 



G(x,y) 



2tt 

y 



ln2 



for x < 1, (31) 



where Q = Q2 ^ the scaling function defined in Eq. (27) 
(x = T//x) [36]. Without the additive constants — 1/8 + 
(ln2)/4 and -1/2 + (ln2)/2, Eqs. (30) coincide with the 
mean-field result assuming an effective interaction con- 
stant g(jn). These constants can be omitted in the weak- 
coupling limit g(jji) C 1. As pointed out in Refs. [9, 10], 
in the weak-coupling limit - where the BKT transition 
temperature Tbkt can be easily determined (see Sec. VI) 
- the approximation (31) remains remarkably accurate all 
the way down to the transition point (/jl/T)bkt- We also 
observe that the limiting behaviors (28,29) and (30,31) 
are very well satisfied not only in the weak-coupling 
limit [9, 10] but also in the strong-coupling limit where 
g(T) > 1 (see Fig. 3). 

The crossover regime |/i| <C T is more difficult to ana- 
lyze in simple terms, and a full numerical solution of the 
RG equations is necessary (see however Sec. IV B and 
Appendix B for an analytical solution in the case /i = 0). 

In the weak-coupling limit g = 2mg <C 1, the scatter- 
ing length <22 is exponentially small. This implies that the 
renormalized interaction constant g(T) ~ g is nearly tem- 
perature independent except for exponentially small tem- 
peratures (which are experimentally unreachable) [ ] . It 
follows that the phase-space pressure and density and the 
entropy per particle, 



V^T)=^(^g), 



(32) 



<S(/i,T) = 2 



V(^T) 



T' 



can be considered as functions of fi/T only, with the 
microscopic interaction constant g entering the scaling 
function J 7 as a parameter. The equation of state then 
exhibits an approximate scale invariance (with no char- 
acteristic energy scales other than /i and T) [9, 10]. 



B. F(0,y)vsy 

The limit \x\ <C 1 is particularly interesting as it cor- 
responds to the quantum critical regime. In this section, 
we discuss the function F(0,y). Figure 4 shows P(0,T), 
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FIG. 4. (Color online) Phase-space pressure V(0, T), phase- 
space density D(0,T) and entropy per particle <S(0,T) vs 
g(T). The symbols show the the data of the ENS, Chicago 
I and Chicago II (with fi = /z c ) experiments [18-20] (see 
Sec. V). 



£>(0,T) and 5(0, T) as a function of g(T). We show in 
Appendix B that 

7>(0, T) ee J-(0, g(T)) - y - ^ In 2 ( -^ ) (33) 

for g(T) -+ 0. The result lim T ^ V(0, T) = tt 2 /6 is exact. 
Experimentally, however, this limiting behavior cannot 
be observed due to the logarithmic temperature depen- 
dence of g(T). In the weak-coupling limit, g(T) = g 
is nearly temperature independent, and the phase-space 
pressure takes the form 



V(0,T)=T(0,g), 



(34) 




FIG. 5. (Color online) (Red) solid lines: dimensionless inter- 
action constant #bh (T) in the two-dimensional Bose-Hubbard 
model for U/t = 0.22 (top) and U/t = 6.25 (bottom) (note the 
different scales on the vertical axes). The (blue) dash-dotted 
lines show the universal limit g(T). 



where ^(O,*?) < lim y _^o ^(0, y) — k 2 /6. In the strong- 
coupling limit, V(0,T) exhibits a weak temperature de- 
pendence coming from that of g(T), but again reaching 
the limiting value \im.T^o'P(0,T) = 7r 2 /6 requires ex- 
tremely small (unrealistic) temperatures. 



C. Thermodynamics of the Bose-Hubbard model 

In this section we discuss the results obtained in the 
two-dimensional Bose-Hubbard model [Eq. (8)] using 
the lattice version of the NPRG [27, 43]. The energy- 
dependent interaction constant is defined by 



#bh0) 



U 



i + c/n(e)' 



with 



n(e) = i 2 1 



d 2 q 



(2tt)2 2(e c 



(35) 



(36) 



where e q is the lattice dispersion of the boson [Eq. (9)]. 
This definition, which is also that used in Ref. [20], 
is justified in Appendix C. In the low-energy limit 
e — »• 0, it coincides with the universal form l~ 2 g{e) 
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FIG. 6. (Color online) Phase-space pressure V(fi c ,T) vs T/t 
for U/t = 6.25. The (green) dashed line shows the result ob- 
tained for a flat density of states (DOS) in the energy window 
[0, St] . The (blue) dash-dotted line shows the universal limit 
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FIG. 7. (Color online) Pressure P(fi c ,T) vs T/t for 
U/t — 6.25. The (green) dashed line shows the result ob- 
tained for a flat density of states in the energy window 
[0, St]. The (blue) dash-dotted line shows the universal limit 

T 2 /(4ntl 2 )T(0rg(T)). 
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FIG. 8. (Color online) Phase-space density D(/i c ,T) ((red) 
solid line) and entropy per particle S(/jl c ,T) ((blue) dash- 
dotted line) vs T/t for U/t = 6.25. The (green) dashed lines 
show the corresponding universal limits. 
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FIG. 9. (Color online) Compressibility k(/jl c ,T) vs T/t for 
U/t — 6.25. The (green) dashed line shows the corresponding 
universal limit. (The numerical noise in the NPRG result 
follows from taking the second-order derivative of the pressure 
with respect to ji.) 



[Eq. (20)], obtained from the continuum model with bo- 
son mass m = l/2tl 2 and scattering length <22 given 
by Eq. (10). The dimensionless interaction constant 
^bh(T) = 2ml 2 g B n(T) = g B n(T)/t is shown in Fig. 5 
for U/t = 0.22 and U/t = 6.25. In both cases, g BB (T) is 
well approximated by its universal limit g(T) [Eq. (20)] 
for T < St (see insets in Fig. 5). 

Figure 6 shows the phase-space pressure V([i c ,T) vs 
T/t for U/t = 6.25 and /j, = \i c = —At. We observe a 
maximum around T/t ~ 2.5 due to the enhanced den- 
sity of states of the square lattice near the band cen- 
ter [ ] . This maximum disappears if we consider a flat 
density of states in the energy window [0, St]. Com- 
paring V([ic,T) and J r (0,^(T)) (with g(T) the universal 
limit of g BB (T) discussed above) we see that the univer- 
sal limit, where V{jjl c -,T) becomes a universal function of 
ma^T, is reached only at very low temperatures T <C t. 
The identification of t as the crossover temperature scale 
for quantum critical behavior is confirmed by the T de- 
pendence of the pressure. For T < t, one finds that 
P(Hc,T) = T 2 /(47rtl 2 )V(fi c ,T) is well approximated by 
the universal limit T 2 /(47rt/ 2 )J r (0,^(T)) (Fig. 7). The 
phase- space density X>(/i c ,T) and entropy per particle 
<S(/i c , T) are shown in Fig. 8 (the low-temperature regime 
where £>(/i c , T) and <S(/i c , T) coincide with their universal 
limits is not shown). 

The fact that the universal regime is reached only 
at low temperatures can also be seen in the tempera- 
ture dependence of the compressibility n = d 2 P/dji 2 
(Fig. 9). Although it is difficult to numerically com- 
pute the second-order derivative of the pressure with re- 
spect to /i, our results clearly show that ft(/i c , T) is below 
the universal limit (1/ Aittl 2 )^ 2 ^ (Q,g(T)). We also note 
that while ft(/i c , T) varies weakly with T in the tempera- 
ture range [£, lOt], it should eventually diverge as T — » 
(see Eq. (B14) in appendix B). We thus disagree with 
the conclusion of Ref. [ ] that quantum criticality is ob- 
served below a characteristic temperature of the order of 
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FIG. 10. (Color online) Phase-space pressure *P(/i, T) vs S/jl/T 
for U/t = 0.22 (top) and U/t = 6.25 (bottom). T/t = 1/100 
(circles), 1/10 (squares), 1 (diamonds) and 3 (triangles). The 
solid lines correspond to F{8^/T,g{T)). For t/U = 0.22, 
fi/T ~ 0.15 at the BKT transition (see table I and Sec. VI). 
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FIG. 11. (Color online) Phase-space pressure V{\l c -> T), phase- 
space density P(/i c ,T), and entropy per particle «S(/i c ,T) vs 
g(T) (U/t — 6.25). The lines show the universal limit ob- 
tained from the scaling function J r (0,^(T)) and its deriva- 
tives. 



the single-particle bandwidth St [47]. 

Figure 10 shows the phase-space pressure V(fJL,T) ver- 
sus 5/i/T for U/t = 0.22 and U/t = 6.25 (Sfi = /i - /i c ), 
and various temperatures ranging from t/100 to 3t. In 
the weak-coupling limit U/t = 0.22, g BB (T) = U/t 
is nearly temperature independent in the temperature 
range [t/100, 3t] (see Fig. 5). At very low temperatures, 
we obtain a perfect agreement between V([i,T) and the 
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universal scaling function F(5ii/T,g) (with g = 0.22). 
At higher temperatures, when T ~ t, we observe that 
V([i,T) slightly deviates from T(8\i/T^ g), in particular 
for large values of \i. This agrees with the previous obser- 
vation that V(nc->T) reaches the universal limit only for 
T < t (see Figs. 6 and 7). For U/t = 6.25, we again find a 
good agreement between V^jl^T) and F(5/j,/T,g(T)) at 
low temperatures and small chemical potential 5/n, but 
deviations are clearly visible at higher temperatures or 
larger values of Sfi/T. 

The phase-space pressure P(/i c ,T), phase-space den- 
sity V(ji c1 T) and entropy per particle <S(/i c ,T) vs g(T) 
are shown in Fig. 11 for U/t = 6.25. The maximum 
around g(T) ~ 3.5 is due to the enhanced density of 
states of the square lattice near the band center [45]. 
For g(T) < 2.5, we recover the universal limit where the 
thermodynamics is determined by the scaling function 
T(8»/T,g(T)). 



V. COMPARISON WITH EXPERIMENTS 

In this section, we compare our theoretical results for 
the scaling function T = T^ with three recent experi- 
ments on two-dimensional Bose gases. The first experi- 
ment was realized with a gas of 87 Rb atoms with scat- 
tering length as = 5.3 nm and a thickness l z = 240 nm 
in the confining direction leading to a dimensionless in- 
teraction constant g = 2mg = 0.22 [19, 48]. The second 
one was performed with 133 Cs atoms and a scattering 
length as controlled by a Feshbach resonance and vary- 
ing in the range 2 — 10 nm resulting in g = 0.1 — 0.52 [18]. 
The last one was realized with a 133 Cs atom gas in 
an optical lattice and can be described by the Bose- 
Hubbard model with t = 2.7 nK, U = 16.7 nK (i.e. 
U/t = 6.25), and a temperature varying in the range 
5.8 - 32 nK (i.e. 2.15* - 32t) [20, 49]. This leads to 
a temperature-dependent dimensionless interaction con- 
stant <?bh(^) varying between 3.95 and 5.75. We refer 
to these experiments as the "ENS", "Chicago I" and 
"Chicago II" experiments, respectively. 

In Fig. 12, we compare the NPRG results with the 
ENS experiment. For g = 0.22, the temperature depen- 
dence of g(T) is negligible so that we expect the scal- 
ing forms (32), which express V, V and S as universal 
functions of ji/T and g, to be very well satisfied. We 
find a nearly perfect agreement between the experimental 
data and the NPRG calculation of the universal function 
T(\i/T^g) (without any fitting parameter). 

In Sec. IV C, we have shown that a Bose gas in an op- 
tical lattice, described by the Bose-Hubbard model with 
U/t = 6.25, reaches the universal limit only at tem- 
peratures of the order of t. In the Chicago II experi- 
ment, the lowest temperature T ~ 2.15t is above t, and 
we should therefore expect experimental data to agree 
only approximately with results obtained from the uni- 
versal function T . The temperature dependence of the 
phase- space density V(ji Cl T) is shown in Fig. 13. There 




FIG. 12. (Color online) Phase-space pressure V(fi,T), phase- 
space density T>(/j, T) and entropy per particle «S(/i, T) vs \i/T 
in the ENS experiment [ ]. The (red) solid lines show the 
NPRG results. 
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FIG. 13. (Color online) Temperature dependence of the 
phase-space density T>(fi c ,T) in the Chicago II experiment 
(triangles) [ )]. The (red) solid line shows the NPRG result 
and the (green) dashed one the universal limit. 
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FIG. 14. (Color online) Phase-space pressure V(/j>, T) vs Sfi/T 
for T = 6.7 nK and T = 11 nK in the Chicago II experi- 
ment [20] . The solid lines show the NPRG results obtained in 
the Bose-Hubbard model and the dashed lines the universal 
limit. 
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FIG. 15. (Color online) Phase-space density P(/x, T) and en- 
tropy per particle S(fi,T) vs Sfi/T for T = 6.7 and T = 11 nK 
in the Chicago II experiment [ ] . The solid and dashed lines 
show the NPRG results obtained in the Bose-Hubbard model. 



is an overall agreement between the experimental data 
and the NPRG results but the existence of a plateau 
for T < St followed by a strong suppression of V(p c ,T) 
at higher temperatures, as advocated in Ref. [20], is 
not supported by the theory. In Fig. 14 we show the 
phase-space pressure V(p,T) vs Sp/T for T/t = 2.5 and 
T/t = 4.1. As expected the NPRG results show devi- 
ations from the universal limit F(5/j,,g(T)) ( note that 
g(T) is nearly temperature independent in the tempera- 
ture range [2. 5t — 4.1t]). For large and negative chemical 
potential 5fi, the pressure is very well approximated by 
the classical dilute gas expression 



7>(/i,T)=4^te-l^l/ T / 2 



L 



,/T 



(37) 



The difference with the universal limit V = e _ ' 5/x l/ T 
(Sec. IV A) is entirely due to the difference between the 
lattice dispersion e q [Eq. (9)] and the free quadratic dis- 
persion q 2 /2ra with m = l/2tl 2 . For T/t = 4.1, Eq. (37) 
gives V ~ 1.3e-l 5/i l/ T when Sp/T < -2. On the other 
hand the experimental data show a remarkable agree- 
ment between the phase space pressure V and the uni- 
versal scaling function T with only a small difference for 
positive dp. Such an agreement is difficult to understand 
in the framework of the Bose-Hubbard model. In par- 
ticular, one would expect V to differ from e~^^l T for 
large and negative Sp and T/t c^ 2 — 4 due to lattice ef- 
fects (see the discussion above). The phase-space density 
V(p, T) and the entropy per particle <S(/i, T) vs Sp/T for 
T/t = 2.5 and T/t = 4.1 are shown in Fig 15; there is a 
good agreement between theory and experiment. 

The ENS, Chicago I and Chicago II experiments can 
be used to obtain P(/x c ,T), V(p c ,T) and <S(/i c ,T) as 
a function of the effective interaction constant g(T). 
The results are shown in Fig. 4. For all three experi- 
ments, we obtain a very good agreement with the uni- 
versal limit (23,24). This confirms that both the ENS 
and Chicago I experiments deal with a weakly inter- 
acting Bose gas in the universal regime. As for the 
Chicago II experiment, such a good agreement is par- 
tially accidental since for #bh(T) — d(T) — 4.3 (the rel- 
evant value of #bh(T) corresponding to the experimen- 
tal data shown in Fig. 4), S(p Cl T) turns out to be very 
close to the universal limit even though the system has 
not reached the universal regime yet (see Fig. 11). We 
also note that for this value of this interaction constant, 
V(p c ,T) and P(/i c ,T) are nearly equal, which implies 
that P(/i c , T) ~ TD(p c1 T) as observed in the Chicago II 
experiment. 



VI. BKT TRANSITION TEMPERATURE 

In this section we show how the BKT transition tem- 
perature Tbkt can be estimated from the NPRG ap- 
proach. For the classical 0(2) model, the NPRG re- 
produces most of the universal properties of the BKT 
transition [50, 51]. In particular one finds a value p$ of 
the dimensionless order parameter (the spin-wave "stiff- 
ness") such that the beta function /3(po,fc) — kdkPo,k 
nearly vanishes for po,/c > Po (here k denotes the RG 
momentum scale, see Appendix A). This implies the ex- 
istence of a line of quasi- fixed points and enables to iden- 
tify a low-temperature phase (T < Tbkt) where the run- 
ning of the stiffness po,k, after a transient regime, be- 
comes very slow, implying a very large (although not 
strictly infinite as expected in the low-temperature phase 
of the BKT transition) correlation length £. In this low- 
temperature phase, the anomalous dimension rjk depends 
on the (slowly varying) stiffness po,/c- It takes its largest 
value ~ 1/4 when the RG flow crosses over to the disor- 
dered (long-distance) regime (for po,/c ~ Pq anc ^ ^ ^ £ _1 )> 
and is then rapidly suppressed as /?o,fc further decreases. 
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FIG. 16. (Color online) Flow trajectories in the plane (n s ,rj) 
for a two-dimensional Bose gas with a dimensionless interac- 
tion constant g — 0.22. The vertical line indicates the value 
of n* . The (red) solid line shows the line of quasi-fixed points 
for n s > n*. The critical trajectory (which joins the line of 
quasi-fixed points for n s = n*) corresponds to fi/T ~ 0.154. 



On the other hand, the beta function is well approxi- 
mated by /3(p 0j fc) = const x (pg - p ,/c) 3/2 for p 0jk < pg, 
and the essential scaling £ ~ e const/(T-T BK T) f ^q 
correlation length above the BKT transition tempera- 
ture Tbkt is reproduced [51]. Thus, although the NPRG 
approach does not yield a low-temperature phase with 
an infinite correlation length, it nevertheless allows us to 
estimate the BKT transition temperature from the value 
of Pq . A reasonable estimate of the BKT transition in the 
two-dimensional XY model has been obtained using the 
lattice NPRG [ ]. Here we use the NPRG to determine 
the BKT transition temperature in a two-dimensional 
Bose gas [53]. 

The flow trajectories in the plane (n 3 ,rj) are shown in 
Fig. 16 for the continuum model with g = 0.22. n s de- 
notes the superfluid density and is analog to the dimen- 
sionless order parameter po of the classical 0(2) model. 
At sufficiently low temperatures, the trajectories join a 
line of quasi-fixed points where the RG flow is very slow, 
before eventually crossing over to the disordered phase 
(n s — >• 0). The value of n s at the merging point with 
the line of quasi- fixed points depends on the temperature 
and chemical potential of the Bose gas. We estimate the 
BKT transition temperature by the trajectory for which 
the merging point corresponds to the value nj (analog 
to Pq i n the classical 0(2) model) of the superfluid den- 
sity. A precise determination of the value of n* (which 
can be obtained by fitting the beta function kdkn Sj k by 
const x (n* — n^) 3 / 2 for n 8i k < ti*) is however difficult as 
it requires the full calO(d 2 ) expansion of the effective ac- 
tion while we solve the NPRG equation within a simple 
truncation of the effective potential [Eq. (Bl)]. Never- 
theless, since the BKT transition in the Bose gas model 
and the classical 0(2) model is controlled by the same 
fixed point, we expect the ratio n^/nf^, where n™ ax is 
the value of n s for which 77 is maximum (Fig. 16), to be 
equal to Pl/p™*. 
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FIG. 17. (Color online) Ratio (/x/T)bkt vs g(T) for T = 
Ta/10 and T = T A /50, where T A = A 2 /2m. The (green) 
dotted line corresponds to the expression (40) with the Monte 
Carlo result C = 13.2 [9, 10]. 



Using this method, we have verified that the ratio p/T 
at the BKT transition is a universal function of g(T), i.e. 



vi /BKT 



H(g{T)), 



(38) 



with % a universal function. Equivalently, since g(T) is 
a function of ma^T, (p/T) bkt can be seen as a universal 
function of ma^T or ma^\i. Figure 17 shows (/x/T)bkt 
obtained for two different temperatures, T = Ta/10 and 
Ta/50 (Ta = A 2 /2m), and a range of values of g. The 
universal form (38) is well satisfied in the weak-coupling 
limit (g(T) ~ g < 1). In this limit, we find 
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0.982 

2tt 



gin 



2 x 9.48 
9 



(39) 



in good agreement with the weak-coupling result [9, 10, 

12, 35] 

Vi /BKT 27T \ 9 J 

where £ c± 13.2 ±0.4 has been obtained from Monte Carlo 
simulation [ ,10]. We ascribe the violation of universality 
at strong coupling, as seen in Fig. 17, to a poor descrip- 
tion of the BKT transition by the NPRG when g>\ [54]. 
Although we can use the same method to determine the 
BKT transition temperature in the Bose-Hubbard model, 
we cannot compare with the experimental result of the 
Chicago II experiment [2 0] which corresponds to a strong- 
interaction regime (<?bh(T) ~ 4.3) where this method is 
not reliable. 



VII. CONCLUSION 

The scale invariance of the equation of state of a 
weakly interacting Bose gas, i.e. the fact that the phase- 
space pressure V{ii,T) depends only on ji/T when the 
dimensionless interaction constant g is small, is well 



13 



understood both experimentally and theoretically. We 
have shown that, more generally, the phase- space pres- 
sure 7^(/i, T) is a universal function of ji/T and the 
temperature-dependent dimensionless interaction con- 
stant g(T) [Eq. (1)]. Using the NPRG approach, we have 
computed the corresponding universal scaling function 
T{x, y) for a two-dimensional gas from weak to strong 
coupling. Recent measurements of the pressure, density 
and entropy in a weakly two-dimensional Bose gas [18, 19] 
allow us to determine both the x and y dependence of 
J~{x, y) in some limits, and the results are found to agree 
remarkably well with the NPRG predictions. 

We have also compared our theoretical results in the 
Bose-Hubbard model with recent experimental data ob- 
tained in a two-dimensional Bose gas in an optical lattice 
near the vacuum-superfluid transition [20]. Our theoreti- 
cal analysis shows that the lowest temperature (T = 2.5£) 
reached in the experiment remains slightly above the 
crossover temperature T ~ t to the quantum critical 
regime where the thermodynamics is fully determined 
by the universal scaling function T . However, some- 
what surprisingly, the experimental data do not show the 
small deviations from (universal) quantum critical behav- 
ior that are expected for T = 2.5£ (see the discussion in 
Sec. V). 

The experiment reported in Ref. [ ] shows that it 
is now possible to measure the thermodynamics of a 
two-dimensional Bose gas in an optical lattice near the 
superfluid-Mott-insulator transition (where the Mott in- 
sulating phase is not the vacuum). Since this transi- 
tion (when it is induced by a density change) belongs 
to the dilute Bose gas universality class, the thermody- 
namics in the superfluid phase is also determined by the 
scaling function T (the BKT transition temperature be- 
ing determined by the scaling function H, see Eq. (38)). 
The nonuniversal parameters m and a<2 should be un- 
derstood as the effective mass and effective scattering 
length of the elementary excitations at the (nontrivial) 
QCP between the superfluid phase and the Mott insula- 
tor. We have recently shown that Eq. (1) indeed holds for 
a three-dimensional Bose gas in an optical lattice near the 
Mott transition and computed the non-universal param- 
eters m and as in the framework of the Bose-Hubbard 
model [55]. Measuring the thermodynamics near the 
superfluid-Mott-insulator transition of a two- or three- 
dimensional Bose gas would allow for a very strong test 
of universality in strongly interacting quantum fluids. 
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Appendix A: Non-perturbative RG 

In this section, we briefly review the NPRG approach 
to interacting boson systems. The strategy of the NPRG 
is to build a family of models indexed by a momentum 
scale k varying from a microscopic scale A down to 0. In 
practice, this is achieved by adding to the (Euclidean) 
action S of the model a term 
rP 
AS k [r^}= dTV>* (q)# fc (qMq). ( A1 ) 

Here V^ 7 ") is a bosonic field and r 6 [0, 1/T] an imagi- 
nary time. The so-called cutoff function iifc(q) vanishes 
for k = so that the action S + ASk=o reduces to the 
action of the model we are interested in. When k is finite, 
i?fc(q) suppresses fluctuations with momenta |q| < k and 
acts as an infrared regulator term. Its value at the micro- 
scopic scale A must be chosen such that the model with 
action S + AS a is exactly (at least numerically) solvable. 
In the standard implementation of the NPRG [56, 57], 
one ensures that all fluctuations are frozen by the ASa 
term, so that the mean-field approximation becomes ex- 
act. For a lattice model, such as the Bose-Hubbard 
model, one can instead choose it^v(q) (with A of the order 
of the inverse lattice spacing) such that t q + Ua(q) — 0. 
The action S + AS a then describes a system of decoupled 
sites (vanishing hopping amplitude) and is exactly solv- 
able. This implementation of the NPRG, referred to as 
the lattice NPRG, was introduced in Ref. [ . } ] and used 
to study the Bose-Hubbard model in Refs. [27, 43]. In 
both the standard and lattice NPRG schemes, the ef- 
fective action of the original system (with action S) is 
deduced from that of the reference system (with action 
S + AS a) by solving a RG equation. 

The main quantity of interest in the NPRG approach 
is the scale-dependent effective action 

T k [^^]= -hiZ k [J*,J]+ f dr [d d r(J*0 + c.c.) 

-AS fc [0*,0], (A2) 

defined as a modified Legendre transform of 
— In Zfc [J* , J] which includes the subtraction of 
ASk[4>*,(f>]- Here J(r,r) is an external (complex) 
source which couples linearly to the boson field ^(r,r), 
Zk [J* , J] the partition function corresponding to the 
action S + AS&, and 



(r,r) 



6 In Z k [J*, J] 
5J*(t,t) 



(V(r,r)) 



(A3) 



the superfluid order parameter. 

In the standard NPRG implementation, the initial 
value T a [</>*,</>] = S[(j)*,4)] of the scale-dependent effec- 
tive action is given by the microscopic action. In the 
lattice implementation (Bose-Hubbard model), 



r A [</>*,</>] = ri oc [</>*,< 



Jo 



dr 



q 



( (q)V(q), (A4) 
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where ri oc is the Legendre transform of the thermody- 
namic potential — In Z\ oc [J* , J] in the local limit (vanish- 
ing hopping amplitude). The effective action T^ can be 
deduced from T\ by (approximately) solving the exact 
flow equation [58] 



dkFkWA] 



.(2) 



> 



d k R 



'fc-fifc 



(r^[<p\4>] + Rk) '}, (A5) 



where T k is the second-order functional derivative of T^. 
In Fourier space, the trace in (A5) involves a sum over 
momenta and frequencies as well as the two components 
of the complex field (f). We refer to Refs. [27, 39, 55] for 
a detailed discussion of the approximations used to solve 
Eq. (A5). 

All thermodynamics properties can be obtained from 
the effective potential defined by 



Vfc(n) 



1 

w 



T k [4>*,4>] 



(A6) 



is a con- 

2 



(V denotes the volume of the system), where 
stant (uniform and time-independent) field and n = \(j)\ 
Its minimum determines the condensate density no,/c and 
the thermodynamic potential (per unit volume) Vo,fc = 
Vk( n o,k) i n the equilibrium state. The pressure is then 
simply given by 



P(H,T) = -V , 



k=0- 



(A7) 



Appendix B: NPRG equations for \i — 



In this section, we discuss the solution of the \i = 
NPRG equation for a continuum model with Hamilto- 
nian (3) (d = 2). We use the cutoff function it^(q) = 
(tk — e q )0(e/c — e q ) (Q(x) denotes the step function) and 
approximate the effective potential by 



Vk(n) = Vb,fc + 5 k n + 



A 



k 2 

n . 



(Bl) 



The initial condition given by the mean-field solution is 
5 a = — [*> and A a = g. The RG equations read 



k = ^A / cot 



Rk(ci) 



2T 



dk^k = 



2 



1 



coth 
2 



T 



sinh 



e q + 4 + i?fe(q" 

+ gfc + flfc(q) 
2T 

2 (eq + fa+Rkfa) 

2T 



(B2) 



where &k = (dkRk)d/dR k and we use the notation J = 
(27r)~ 2 f cPq. Note that with the truncation (Bl), the 
NPRG equations in the normal phase (vanishing con- 
densate density) reduce to standard one-loop equations. 



Except for minor differences (due to a slightly different 
RG scheme), Eqs. (B2) are equivalent to the one-loop RG 
equations derived in Refs. [12, 59] using a momentum- 
shell one-loop RG approach with a sharp cutoff. 

Since e q + i?fc(q) > e&, the RG equations can be ap- 
proximated by their T = limit when e/~ ^> T, 



dkh = 0, 
dk^k 






(B3) 



There is a quantum-classical crossover when k becomes of 
the order of fc^ = v2mT. In the classical regime where 
k <C &t> e q + 5k + Rk(q) <C T for the values |q| ~ k 
contributing to the momentum integrals in (B2), we find 



dkSk = 
dk^k — 



2mT , 

7T 

20m 2 T 

nk 2 



-A* 



(B4) 



^/c' 



To obtain an approximate solution of the RG equations, 
we first integrate (B3) between k = A and k = &t, which 



gives 



4 T ^ 0, 
A, T ~#(T), 



(B5) 



where g(T) = g(T)/2m is defined by (20). We then in- 
tegrate Eqs. (B4) between fc T and k = \/2m5k (the RG 
flow stops beyond this point) with boundary values at kr 
given by (B5). We deduce the approximate expressions 
of 5 = 5 k =o and A = A fe = , 



5 
A 



r ^ln(fl(A + 2mBS) 

2m5 



(B6) 



A + 2ra££ ' 



where 



10ra 2 T 



5 



A 



(/(T) 2mT 



(B7) 



The one-loop RG equations are essentially exact in the 
quantum regime k ^> kr (which coincides with the vac- 
uum limit when /i = 0) but requires A& = 2m\k < Xk T ^ 
1, i.e. g(T) <C 1, to be valid in the classical regime. Using 
k<T and anticipating that 5/T ^> g(T), the equation 
for 5 simplifies into 



5_ 
T 



2tt \T 



The solution of Eq. (B8), 

S = 9(T) 
T 2tt 



"&)• 



(B8) 
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FIG. 18. (Color online) 5/T vs g(T) obtained from the numer- 
ical solution of Eqs. (B2) with /x = ((red) solid line). The 
(blue) dash-dotted line shows the approximate solution (B9) 
with an overall multiplicative constant a — 0.87. The sym- 
bols show 5/T in the Bose-Hubbard model for U/t = 3.33 and 
10 (and ii — fi c ). 



can be written in terms of the Lambert function W(x) 
defined by We w = x. Using W(x) ~ Inx for x ^> 1, we 
obtain 



8_ = m j2 1 _\ 
T 27T \g(T)) 



(B10) 



for g(T) — >> 0. The gap 5 = 5^=o obtained from the 
NPRG equations is shown in Fig. 18. There is a very good 
agreement with the expression (B9) up to a multiplica- 
tive constant a ~ 0.87 (which accounts for the rather 
crude treatment of the quantum-classical crossover when 
solving the RG equations). Figure 18 also shows the gap 
computed in the Bose-Hubbard model. At low temper- 
atures, we recover the universal limit described by the 
continuum model. We also observe a non-monotonous 
variation of 5/T which is due to the enhanced density 
of states of the square lattice near the band center [45]. 
The maximum of 5/T with respect to g(T) is responsible 
for the maximum observed in the phase-space pressure 
V(/jL c ,T) as a function of g(T) (see Fig. 6). 

To obtain the pressure, we adopt the standard 
momentum-shell RG point of view where — <5 is inter- 
preted as a renormalized chemical potential. Comput- 
ing the pressure directly from the renormalized parame- 
ters [59], and ignoring the renormalized interaction A <C 
1, we obtain 



V(0,T)=U 2 (e- 5 / T ), 



(Bll) 



where Li2 is a polylogarithm function. Together with the 
condition 5 <C T and the asymptotic behavior of Li2(x) 
for \x - 1| < 1 and x < 1, Eq. (Bll) gives 



V(0,T) 



7T 



^ \9(T)J 



(B12) 



This expression agrees with the numerical solution of 
Eqs. (B2) but differs from the analytical result reported 



in Ref. [ ]. The result lim T ^ ^(0,T) = tt 2 /6 is exact 
since the one-loop approximation becomes exact in the 
limit g(T) -► 0. 

To compute the density D and the compressibility ft, 
we start the RG procedure with an infinitesimal chemi- 
cal potential, i.e. 5\ = —ll. Integrating the RG equa- 
tions in the classical and quantum regimes, we then find 
5k T — —[i and 5 (/a) = 5(0) + 5k T , where —5(0) is the 
renormalized chemical potential for \i = [Eq. (B9)]. 
We deduce 



D(0,T) 



k(0,T) 



dP(j*,T) 



d 2 P(v,T) 



n=o 



dfi 2 



fj,=0 



dP(0,T) 

OS ' 

<9 2 P(0,T) 
W 2 ' 



(B13) 



Using 1a' 2 {x) = — ln(l — x)/x, we finally obtain 



p(0,T) = -ln(l 



«(0,T) 



m 



1 m 



(B14) 



2n e s ' T - 1 



5(T)ln(^)' 



a result which is also in good agreement with the numer- 
ical solution of the RG equations. Note that the com- 
pressibility ft(0, T) diverges at the quantum critical point 
/i = T = 0, in agreement with the known result 



m 

2^ 



fr(/i, T = 0) ~ — £- ln ( - ^jvna\\i 



2m 

#0) 



(B15) 



for small positive \i (see Eq. (26)). 



Appendix C: Interaction constant gBii(e) in the 
Bose-Hubbard model 



In the continuum model, the energy-dependent in- 
teraction constant is defined from the solution of the 
T = Li — RG equation. The latter can be written 

as 



9(e) = 



1+ 511(e) ' 



(CI) 



where 



duo 

2^ 



©(^q - e) 



J q J — ( 

= /e(e q -e)-L 

^q Zfc q 



[iuj - e q )(-iw - e q ) 



(C2) 



is a particle-particle propagator with an infrared cutoff 
e. Equation (CI) makes it clear that the one- loop RG 
equation resums the ladder diagrams contributing to the 
two-particle vertex. For bosons with a quadratic disper- 
sion e q = q 2 /2m (and an energy-independent density of 
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states N(e) = f S(e — e q ) = m/2ir), Eq. (C2) is equiva- 



lent to 



n(e) = / q ^ 



+ e) 



(C3) 



in the limit e <C A 2 /2m. In the Bose-Hubbard model, we 
define the energy-dependent interaction constant #bh(<0 
from (CI) and (C3) with g = U and e q the lattice disper- 
sion (9) [see Eqs. (35,36)]. Note that using (C2) rather 
than (C3) would yield the same universal limit (20) in 
the low-energy limit. 



[1] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957); 
T. D. Lee, K. Huang, and C. N. Yang, Phys. Rev. 106, 
1135 (1957). 
[2] S. Sachdev, Quantum Phase Transitions (Cambridge 

University, Cambridge, England, 1999). 
[3] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
[4] Q. Zhou and T.-L. Ho, Phys. Rev. Lett. 105, 245702 

(2010). 
[5] K. R. A. Hazzard and E. J. Mueller, Phys. Rev. A 84, 

013604 (2011). 
[6] Although irrelevant in the renormalizat ion- group sense, 
the boson-boson interaction constant necessary enters the 
equation of state since an ideal Bose gas is unstable for 
[i > 0. Furthermore, g(T) is a dougerously irrelevant vari- 
able above the upper critical dimension d^ = 2, which 
means that the scaling function Td [Eq. (1)] may depend 
on g(T) in a singular way. 
[7] E. Braaten and H.-W. Hammer, Physics Reports 428, 

259 (2006). 
[8] In a three-dimensional Bose gas, I is given by the van der 
Waals length / v dw , which is typically of the order of the 
s-wave scattering length 03. In a quasi- two-dimensional 
gas, it is given by the thickness l z of the gas in the con- 
fining direction (see Sec. II). 
[9] N. Prokof'ev, O. Ruebenacker, and B. Svistunov, Phys. 
Rev. Lett. 87, 270402 (2001). 

[10] N. Prokof'ev and B. Svistunov, Phys. Rev. A 66, 043608 
(2002). 

[11] V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1970); 34, 
610 (1971); J. M. Kosterlitz and D. J. Thouless, J. of 
Phys. C 6, 1181 (1973); J. Phys. C 7, 1046 (1974). 

[12] D. S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 
(1988). 

[13] C.-H. Cheng and S.-K. Yip, Phys. Rev. B 75, 014526 
(2007). 

[14] T.-L. Ho and Q. Zhou, Nat. Phys. 6, 131 (2009). 

[15] Z. Hadzibabic, P. Kriiger, M. Cheneau, B. Battelier, and 
J. Dalibard, Nature 441, 1118 (2006). 

[16] P. Clade, C. Ryu, A. Ramanathan, K. Helmerson, and 
W. D. Phillips, Phys. Rev. Lett. 102, 170401 (2009). 

[17] S. P. Rath, T. Yefsah, K. J. Giinter, M. Cheneau, R. Des- 
buquois, M. Holzmann, W. Krauth, and J. Dalibard, 
Phys. Rev. A 82, 013609 (2010). 

[18] C.-L. Hung, X. Zhang, N. Gemelke, and C. Chin, Nature 
470, 236 (2011). 

[19] T. Yefsah, R. Desbuquois, L. Chomaz, K. J. Giinter, and 
J. Dalibard, Phys. Rev. Lett. 107, 130401 (2011). 

[20] X. Zhang, C.-L. Hung, S.-K. Tung, and C. Chin, Science 
335, 1070 (2012). 

[21] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 



[22] D. S. Petrov, M. Holzmann, and G. V. Shlyapnikov, 
Phys. Rev. Lett. 84, 2551 (2000). 

[23] L.-K. Lim, C. M. Smith, and H. T. C. Stoof, Phys. Rev. 
A 78, 013634 (2008). 

[24] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

[25] For a calculation of the scattering length in the two- 
dimensional Bose-Hubbard model, see Ref. p ]. 

[26] Since the QCP ground state is the vacuum, a particle 
added to the ground state propagates as a free particle. 
The single-particle propagagor is therefore given exactly 
by the non-interacting propagator [Eq. (11)]. 

[27] A. Rancon and N. Dupuis, Phys. Rev. B 84, 174513 
(2011). 

[28] See, for instance, S. K. Ma, Modern Theory of Criti- 
cal Phenomena (Advanced Books Classics, New Edition, 
2000). 

[29] The dimensionless variables T, p, and P are defined ac- 
cording to their scaling dimensions at the fixed point: 
[f] = z, \p\ = l/i/ and [P]=d+ z. With the additional 
2ra factors (and the fact that z — \/v — 2 at the Gaus- 
sian fixed point), these quantities have no physical units. 

[30] In the RG terminology, these higher-order interaction 
contants give rise to corrections to scaling. 

[31] When these higher-order contributions are taken into ac- 
count, the scaling function T must involve additional ar- 
guments, such as the dimensionless three-body interac- 
tion constant, and the low-energy physics cannot be de- 
scribed in terms of the mass m and the scattering length 
ad only. In two dimensions however, it turns out that all 
n-body interaction constants (n > 2) can be expressed as 
a function of a 2 (see Ref. [ 12]) and the universal scaling 
function Ti remains a function of ji/T and a<i . 

[32] S. K. Adhikari, T. Frederico, and I. D. Goldman, Phys. 
Rev. Lett. 74, 487 (1995). 

[33] M. Schick, Phys. Rev. A 3, 1067 (1971). 

[34] V. N. Popov, Theor. and Math. Phys. (Sov.) 11, 236 
(1972). 

[35] V. N. Popov, Functional Integrals in Quantum Field The- 
ory and Statistical Physics (Reidel, Dordrecht, Holland, 
1983). 

[36] Noting that g(e) is a function of ma d e, the universal scal- 
ing function Qa can be easily related to Td- 

[37] J. O. Andersen and M. Strickland, Phys. Rev. A 60, 1442 
(1999). 

[38] C. Wetterich, Phys. Rev. B 77, 064504 (2008). 

[39] N. Dupuis and K. Sengupta, Europhys. Lett. 80, 50007 
(2007); N. Dupuis, Phys. Rev. Lett. 102, 190401 (2009); 
Phys. Rev. A 80, 043627 (2009). 

[40] S. Floerchinger and C. Wetterich, Phys. Rev. A 77, 
053603 (2008); Phys. Rev. A 79, 063602 (2009). 

[41] S. Floerchinger and C. Wetterich, Phys. Rev. A 79, 
013601 (2009). 



17 



[42] A. Sinner, N. Hasselmann, and P. Kopietz, Phys. Rev. 

Lett. 102, 120601 (2009); Phys. Rev. A 82, 063632 

(2010); C. Eichler, N. Hasselmann, and P. Kopietz, 

Phys. Rev. E 80, 051129 (2009). 
[43] A. Rancon and N. Dupuis, Phys. Rev. B 83, 172501 

(2011). 
[44] The temperature for which g(T) = g/2, T = 

(2/ma|)exp(-2C- I6ir/g) = (A 2 /2ra)(Aa 2 /2) 2 e 2C , is 

proportionnal to a 2 . 
[45] The density of states of the square lattice is given by 

*«?) = f*8{e-e«) = (2nH)- 1 K[^{^W/4t] (0 < 

e < St) [with e q defined in Eq. (9)] . The complete elliptic 

integral of the first kind K(x) diverges logarithmically for 

x^0. 
[46] S. Fang, C.-M. Chung, P. N. Ma, P. Chen, and D.-W. 

Wang, Phys. Rev. A 83, 031605 (2011). 
[47] Note that the divergence of the compressibility k at the 

quantum critical point jjl = T = was overlooked in 

Ref. [ ]. 
[48] Our definition of the dimensionless interaction constant 

g = 2mg differs by a factor of 2 from that used in 

Refs. [18-20]. 



[49] Because of a nonvanishing next-nearest neighbor hop- 
ping in the Chicago II experiment [ ] , the effective mass 
m slightly differs from l/2£/ 2 (X. Zhang and C. Chin, 
private communication). This difference leads to minor 
modifications of our results and has been ignored. 

[50] M. Grater and C. Wetterich, Phys. Rev. Lett. 75, 378 
(1995). 

[51] G. V. Gersdorff and C. Wetterich, Phys. Rev. B 64, 
054513 (2001). 

[52] T. Machado and N. Dupuis, Phys. Rev. E 82, 041128 
(2010). 

[53] See also Ref. [41]. 

[54] While the beta function /3(n s ,k) = kdkn s ,k for n s ,fc > 
n* is very strongly suppressed at low temperatures in 
the weak- interaction regime, when g(T) > 1 its value 
significantly increases and the NPRG description of the 
BKT transition is not reliable any more. 

[55] A. Rancon and N. Dupuis, Phys. Rev. A 85, 011602(R) 
(2012). 

[56] J. Berges, N. Tetradis, and C. Wetterich, Phys. Rep. 
363, 223 (2002). 

[57] B. Delamotte, cond-mat/0702365. 

[58] C. Wetterich, Phys. Lett. B 301, 90 (1993). 

[59] S. Sachdev, T. Senthil, and R. Shankar, Phys. Rev. B 
50, 258 (1994). 



